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ABSTRACT 

This annual report details progress in basic research in thermoacoustic heat transport made during 
the period October 1, 1989 through September 30, 1990. Four projects are discussed. The first 
project involves measurements of thermoacoustic phenomena using thermoacoustic couples 
(TACs). The purpose of the second project is to examine means of visualizing thermoacoustic 
processes. The third project is an investigation of heat driven prime movers. The fourth project 
considers the performance of thermoacoustic heat engines that operate with an acoustic sound field 
that is not a perfect standing wave. A publications, patents, presentations, and honors report is 
also included. 
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A. INTRODUCTION 



The purpose of this research project is to investigate thermoacoustic heat transport 

phenomena. The premise of the majority of this research can be summarized as follows. While 

early measurements of thermoacoustic phenomena qualitatively agree quite well with the theory 

developed by Wheatley 1 and Swift, there have been few thorough, quantitative comparisons of 

measurement and theory. The quantitative studies that have been performed previous to our 

3 

involvement indicate that the theory does not adequately explain the results. 

Thermoacoustic phenomena and the associated devices can be described in terms of two 
basic modes of operation, heat pumps and prime movers. In thermoacoustic heat pumps, an 
acoustic field stimulates the transport of entropy over the surface of an object, resulting in a thermal 
gradient being established across that object. In prime movers, a thermal gradient is established 
across an object which is part of, or housed in, an acoustic resonator. At sufficiently high 
gradients, sound is spontaneously generated in the resonator. The stored thermal energy is 
converted into acoustic energy; the work output of the device is equal to the energy in the acoustic 
field. Our study of thermoacoustic processes involves investigating both of these fundamental 
classes of phenomena. Our previous research has indicated that when the acoustic amplitude is 
low, measurements and theory agree reasonably well. However, when the amplitude is large, 
there is considerable disagreement. 4 * 8 In other words, in the linear acoustic regime, 
thermoacoustics seems to be relatively well understood, whereas in the nonlinear regime, it is not. 
The question now arises "Is it the thermoacoustics which in not well understood, or the nonlinear 
acoustics?" 

In FY 1990, we extended our previous measurements and began to investigate ways to 
separate apparent inadequacies of the theory of thermoacoustics from phenomena caused by 
nonlinear acoustic affects. In addition, we have begun to consider the performance of 
thermoacoustic heat engines that operate with an acoustic sound field that is not a perfect standing 
wave. Pertinent previous results of our investigations, along with a brief discussion of the 
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progress made during FY 90, follow. 



B. SUMMARY OF PROGRAM 

Our research program consists of four projects. The first project involves measurements of 
thermoacoustic phenomena using thermoacoustic couples (TACs). The purpose of the second 
project is to examine means of visualizing thermoacoustic processes. The third project is an 
investigation of heat driven prime movers. The fourth project considers the performance of 
thermoacoustic heat engines that operate with an acoustic sound field that is not a perfect standing 
wave. Summaries of progress made in these projects during FY 1990 follow. 

1. TAC Measurements 

We have devoted a large effort to a thorough investigation of basic thermoacoustic 
phenomena using the simplest class of thermoacoustic engine - a theimoacoustic couple, or TAC. 
A typical TAC is composed of a stack of short, poorly thermally conducting plates, spaced by at 
least a few thermal penetration depths. The advantage of using a TAC is that its design allows 
assumptions to be made which reduce the theory of thermoacoustic heat transport to its simplest 
form. Our previous research 4- * 5 has shown that the basic theory agrees relatively well with the 
result of measurements made at low acoustic intensities. However, there is serious disagreement 
between theory and experiment at higher acoustic intensities. A major effort prior to FY 90 has 
been to develop a computer controlled apparatus which would facilitate data acquisition. An 
extensive set of measurements has been made with this apparatus, which allows us to vary the 
acoustic pressure amplitude, the position of the TAC, the acoustic frequency, and the mean gas 
pressure and temperature. These measurements represent one of the most complete studies of 
thermoacoustic phenomena to date. 
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The first part of FY 90 was devoted to analyzing and communicating the results of the initial, 
exploratory, phase of the TAC measurements. 9 Several conclusion may be drawn from this work. 
The first conclusion is that the thermoacoustic effect is relatively, yet not completely, well 
understood at low values of the drive ratio (the ratio of the peak acoustic pressure amplitude at a 
pressure antinode to the mean pressure of the gas) for the plate configurations used in these 
measurements. The second conclusion is that there are two distinct regions of behavior at higher 
drive ratios. One region, which occurs for drive ratios less that approximately 1.1%, is 
characterized by a linear decrease in the agreement between theory and experiment. (See Reference 
9 for the meaning of "linear decrease in the agreement.") The second region, which occurs for 
drive ratios greater than approximately 1.1%, is characterized by the onset of irregularities in the 
temperature difference data series. Also, the linear dependence of the disagreement disappears. 
Further, because the agreement is, generally, best in the vicinity of the pressure antinode (velocity 
node), some velocity dependent effect may be the cause of the poor agreement at higher drive 
ratios, where the acoustic particle displacement amplitude can be as much as a few millimeters. 
There is also evidence of an effect which depends on the produce of acoustic pressure and velocity. 

In light of these results, we have initiated the second phase of the TAC measurements, 
intended to further probe the discrepancies and behaviors uncovered during the first phase. The 
first series of measurements is designed to investigate whether or not edge effects become 
important at the higher drive ratios and are the cause of the irregularities discussed above. In the 
previous measurements, we measured the temperature difference developed across the TAC with a 
series of thermocouples (called a thermopile) placed along the front and back edges of the TAC, as 
illustrated in Fig. 1. Because of its proximity to the edge, the thermopile may be sensitive to 
effects which, though perhaps causing local deviations in the temperature profile, do not affect the 
temperature profile in interior regions of the plate. In other words, although we observe 
considerable deviations between theory and measurements made at the plate edges, the interior (and 
majority) of the plate may behave in accordance with predictions. 

It is useful to examine some data taken with a TAC having an edge-to-edge thermopile. An 
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Figure 1: The general design of a thermoacoustic couple (TAC) used in the measurements. The 
central plate is instrumented with a thermopile used to measure the temperature 
difference developed across the TAC. Four guard plates surround the central plate. 




0.5 1 1.5 2 2.5 3 3.5 4 4.5 



kx 

Figure 2: Graph showing the measured (dashed curve) and the theoretical value (solid curve) of 
the temperature difference developed across a TAC as functions of kx for helium and a 
drive ratio of 1.99%. The mean gas pressure is 114.1 kPa. The acoustic pressure 
amplitude is 2.27 kPa and the frequency is 696 Hz. The theoretical temperature 
difference is computed using Eq. (1) of Ref. 9. 
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example is shown in Fig. 2, which illustrates the comparison between the measured temperature 
difference across a TAC (dashed curve) and the theoretical value (solid curve). (Refer to Ref. 9 for 
full details of the TAC construction and experimental parameters.) The TAC is in helium at 298.4 
K. The mean gas pressure and temperature are 114.1 kPa and the acoustic pressure amplitude is 
2.27 kPa, giving a drive ratio of 1.99%. The frequency is 696 Hz. The theoretical temperature 
difference is computed using Eq. (1) of Ref. 9. The temperature differences are plotted as 
functions of kx, where k is the wave number and x is the location of the center of the TAC. The 
positions of the particle velocity and pressure antinodes are indicated on the graph. Four areas for 
comparisons are obvious. First, the agreement is quite good in the vicinity of the pressure 
antinode. The slopes of the measured and predicted temperature differences are approximately 
equal. Secondly, the agreement in the vicinity of the velocity antinode is poor. The magnitude of 
the slope of the measured temperature difference is much less than that of the predicted result. 
Third, the data series of the measured temperature difference is quite irregular and jagged in the 
vicinity of the maximum theoretical temperature difference. Finally, the measured maximum 
temperature difference is less than and occurs at a different position than the predicted maximum. 
In short, the agreement between theory and experiment is poor except near the pressure antinode. 

Figure 3 illustrates the theoretical temperature difference as a function of kx for drive ratios 
from 0.17 to 1.99% in helium at a mean pressure of 1 13.8 ± 0.2 kPa, temperature 298.4 ±0.1 K, 
and frequency 702.8 ± 4.5 Hz. (The individual values of the drive ratio, mean pressure, 
temperature and frequency have been omitted, in order to reduce clutter on the graph. The values 
given above represent the average and standard deviation for the entire data set The data presented 
in Fig. 2 is a subset of these data.) This graph clearly indicates the progression from a sinusoid to 
a sawtooth curve and the displacement of the maximum temperature difference toward the pressure 
antinode as described by Wheatley gt al. 1 The measured temperature differences are shown in Fig. 
4 for the same experimental conditions. The measurements do confirm qualitatively the sinusoid- 
to-sawtooth progression. One of the curves is dashed to demarcate two regions of behavior. The 
curves between the kx axis and the dashed curve, corresponding to lower drive ratios, are smooth 
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Figure 3: Graph showing the theoretical temperature difference as a function of kx for drive ratios 
from 0.17 to 1.99% for a TAC in helium at approximately the same mean pressure and 
frequencies as those in Fig.2. The data presented in Fig. 2 is a subset of these data. 




0.5 1 1.5 2 2.5 3 3.5 4 4.5 



kx 

Figure 4: Graph showing the measured temperature difference as a function of kx for drive ratios 
from 0.17 to 1.99% for a TAC in helium at approximately the same mean pressure and 
frequencies as those in Fig.2. The data presented in Fig. 2 is a subset of these data. 
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and regular. The curves on the other side of the dashed curve, corresponding to higher drive 
ratios, become jagged and irregular. This sudden transition occurs in all of the data sets we have 
taken with TACs having a single thermopile spanning the TAC edge-to-edge. The dashed curve 
corresponds to a drive ratio of 1.03%. 

In order to investigate whether edge effects are the cause of any of these discrepancies, we 
have constructed a TAC with two thermopiles, whose junctions do not lie along the edge, and 
repeated some of the previous measurements. Preliminary results are shown in Figs. 5 and 6, 
which show the temperature difference measured by the two thermopiles as a function of kx. The 
TAC is 1.495 cm long. The two thermopiles are 1.373 cm and 0.682 cm long and centered on the 
TAC. Therefore, the thermocouple junctions lie approximately 0.061 cm and 0.407 cm from the 
edges, respectively. Figure 5 shows the data from the longer thermopile, Fig. 6 from the shorter. 
The TAC is in helium at a mean gas pressure and temperature of 200 ± 1 kPa and 295.9 ±0.1 K. 
The frequency range for the data is 662.8 ± 4.5 Hz. The drive ratio range is 0.20% - 2.70%. The 
dashed curve in both figures corresponds to a drive ratio 0.97%. (Recall that the dashed curve in 
Fig. 4 corresponds to a drive ratio of 1.03%.) Comparing these results with those shown in Fig. 
4, they are strikingly similar. (No direct comparison is possible since the experimental conditions 
are different.) Moreover, the transition to an irregular behavior starts at a 1% drive ratio, as it did 
with the edge-to-edge design. These data show that the temperature profile of the interior of the 
TAC behaves the same as that measured at the edge. Therefore, the irregularities previously 
observed are not an artifact of placing the thermopile along the edge. The irregularities may still be 
the result of the edge (e.g., turbulence generated at the edge), but they are not isolated at the edge. 
They extend at least half way to the center of the TAC. 

We will continue this series of measurements in FY 1991 to further identify the extent of 
these irregularities and determine their cause. We also intend to extend our measurements to 
geometries other than TACs, in order to see to what extent the observed behavior is present in 
more useful thermoacoustic engine stacks. 
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kx 



Figure 5: Graph showing the temperature difference developed across a TAG as measured with 
the longer of the two thermopiles. The TAC is in helium at a mean gas pressure and 
temperature of 200 ± 1 kPa and 295.9 ±0.1 K. The frequency range for the data is 
662.8 ± 4.5 Hz. The drive ratio range is 0.20% - 2.70%. The dashed curve in both 
figures corresponds to a drive ratio 0.97%. 
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Figure 6: Graph showing the temperature difference developed across a TAC as measured with 
the shorter of the two thermopiles. The experimental conditions are the same as those 
in Fig. 5. The dashed curve in both figures corresponds to a drive ratio 0 . 97 %. 




9 



2. Visualization of Thermoacoustics 



It is assumed in the theory applied to our TAC measurements that all parts of the plate behave 
in the same manner, that edge effects and other complicating factors do not exist. Therefore, if 
such factors do occur it is important to know about them. One way to investigate whether or not 
these complications play any significant role in our experiments is to make measurements such as 
those described above. However, these measurements with thermopiles give only the average 
temperature difference across portions of the plate. They do not give the actual temperature 
distribution. Perhaps a better way is to find a method to "visualize" the temperature distribution. 
This notion involves a completely different measurement technique. 

Two methods of visualizing thermoacoustics are apparent. One technique, proposed by us 
previously, is .to utilize a thermal imaging camera in order to photograph, in the infrared, a plate or 
stack of plates exposed to an acoustic standing wave. The result would be a real time image of the 
evolution of the temperature gradients in the plate and could yield information about thermoacoustic 
heat transport which can not be gained using discrete thermocouples. A thermal imaging camera is 
available in the Electro-Optics laboratory in the NPS Physics Department. However, there are 
several disadvantages associated with its use. Because the resolution of these disadvantages could 
prove expensive, we have decided to initially investigate a second, though less sophisticated, 
technique to determine whether visualization is fruitful at all. 

In this second technique, strips of liquid crystal, temperature sensitive, sheets replace, for 
instance, the plates of a TAC. When these strips are exposed to a large amplitude acoustic standing 
wave, a visible color gradient is setup in the strip. We have constructed a transparent, Plexiglas 
resonator in which to conduct the experiment so that the strips are readily visible. The strips have a 
5 °C temperature range over which the color changes from brown to green as the temperature 
increases. When a standing wave is setup, the color change occurs quickly. Within a few 
seconds, the edge of the strip closest the nearest pressure antinode turns green, while the other 
edge turns brown. As time progresses, the green and brown regions expand toward the center of 
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the strip, forming a uniform temperature gradient across the strip. Eventually, the mean 
temperature of the strip rises and the overall color pattern shifts towards green. Finally, the entire 
strip indicates that the temperatures are out of the temperature range of the strip. 

Several important results arise from these preliminary experiments. The extreme edges 
change color first, demonstrating that the transported heat accumulates at the edges of the plates 
first, as first pointed out by Wheatley £i ai. 1 A uniform gradient is established only after thermal 
conduction acts to destroy this initial temperature gradient. These observations are entirely 
consistent with the present understanding of thermoacoustics. A final, yet not unexpected, result is 
that mean temperature of the plate increases. 

Although, we have made only preliminary investigations with the liquid crystal method of 
visualizing thermoacoustics, the method has already yielded useful qualitative information. Aside 
from these benefits, this technique is very useful for demonstrating thermoacoustics. 10 

3. Heat Driven Prime Movers 



The largest effort in FY 90 has been devoted to the investigation of heat driven prime 
movers. There are two distinct regions of operation of a prime mover - below and above the onset 
of self-oscillation. We are investigating each. In an effort to establish a thermoacoustics research 
program at the University of Mississippi, Professor Henry Bass spent his Sabbatical leave at NPS 
to learn from our experience. He has been actively involved in the prime mover research. This 
exchange has not been unidirectional. Dr. Bass's expertise in propagation of sound in porous 
media has been invaluable to our efforts. 

In order to investigate the work output of a prime mover below onset, use is made of the fact 
that the conversion of thermal energy to acoustic energy affects the net absorption of acoustic 
energy in the prime mover. As the temperature difference across the prime mover stack (See Fig. 
7.) increases from zero to its value at onset, the net absorption decreases to zero. The absorption 
coefficient can be determined by measuring the quality factor Q of the prime mover, which in turn, 
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Figure 7: Diagram of a thermoacousric prime mover. 




Figure 8: Graph of 1/Q vs the temperature difference across the prime mover stack filled with 
helium at a mean pressure of 238 kPa. The curve is the results of calculations 
discussed in Ref. 11. 



12 



can be related to the acoustic power dissipated in the tube. The dissipated power can be calculated 
using equations derived by Swift, 2 which express the acoustic power output of a thermoacoustic 
engine as a function of the applied temperature difference across the stack. The absorption 
coefficient can also be extracted from the propagation constant, as is commonly done in porous 
media calculations. 11 

An example of data taken below the onset of self-oscillation is shown in Fig. 8, which is a 
graph of 1/Q, proportional to the net absorption coefficient, as a function of the temperature 
difference across the prime mover stack. 12 The data were obtained by driving the prime mover in 
the vicinity of its fundamental resonance and measuring the frequency response. It is important to 
point out that the acoustic levels involved are well within the linear regime.The curve is based on 
calculations discussed in Ref. 11. Though the overall agreement is quite good, there are some 
discrepancies. There is a tendance to over predict the absorption at low temperature differences. 
Although not so obvious from Fig. 8, other data sets indicate some difficulty in predicting the 
onset temperature The exact reasons for these discrepancies are not entirely clear. What is clear is 
that certain aspects of the theory have not been adequately addressed. 

We have also studied prime movers above the onset of self-oscillation. 13 The prime mover is 
said to have reached "onset" when the temperature difference across the stack is sufficient for the 
prime mover to generate and sustain detectable levels of sound. The use of the word "detectable" 
is not meant to imply that onset is a subtle question of detection thresholds. When onset is 
reached, the observer (and everyone else in the room) knows it. 

Once onset of self-oscillation is reached, the acoustic amplitude in the tube immediately 
assumes a large value, typically about 1% of the ambient pressure. The observed waveform is 
noticeably non-sinusoidal. As more energy is supplied to the hot end of the stack, the temperature 
of that end increases only slightly while the acoustic amplitude in the tube increases rapidly, 
reaching as much as 7% of ambient pressure. Unfortunately, as the acoustic amplitude increases, 
an increasing fraction of the acoustic energy appears as higher harmonics - harmonic distortion 
increases. It is this distortion that has occupied our interest. In order to understand this 
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phenomena, we must not only fully understand harmonic generation from finite amplitude standing 
waves but also any contributing nonlinear generation from the thermoacoustics. 

Figures 9 and 10 show the waveform and the spectrum of the sound generated by the prime 
mover above onset The gas is helium at a mean pressure of 307 kPa. The temperature difference 
across the stack is 325 °C, which is slightly above onset. The signal exhibits slight distortion, 
particularly in the positive cycle. Figure 10 shows that the difference in spectrum level between the 
first few modes is larger than 15 dB. Figures 1 1 and 12 show results for a temperature difference 
of 368 °C, the signal is distorted sharply in both positive and negative half-cycles. The difference 
between spectrum levels for the first few modes has decreased further to less than 6 dB and more 
energy has been spread to higher modes. 

There can be little doubt that the signals at the higher frequencies are being generated by the 
fundamental. If we examine the frequency of these spectral peaks, they are integral multiples of 
the fundamental to within four significant digits and have no measurable line width. Prediction of 
harmonic generation (one method of describing the non-sinusoidal waveform) is not straight 
forward even with a good understanding of dissipation mechanisms in the stack. In a resonant 
tube, the frequencies of the overtones depend upon dispersion while harmonic generation results in 
exact multiples of the fundamental. This means that sound generated nonlinearly will be generated 
at a frequencies different from tube resonances. 

We have applied a treatment of finite amplitude standing waves by Coppens and Sanders 14 
which allows us to calculate the the strengths of the harmonics relative to that of the fundamental. 
Performing the calculations for the situation represented in Figs. 9 and 10 yields (P2/P1) = *20 dB 
and (P3/P1) = -32 dB. Referring to Fig. 10 shows that the measured values are closer to -16 and - 
32 dB, respectively, in reasonable agreement with theory. For the case shown in Fig. 12 the 
calculations predict (P2/P1) = -7 dB and (P3/P1) = -8 dB while the measured values are -7.6 and - 
14.5 dB, respectively. At even higher fundamental amplitudes, the results are less acceptable. 
Although the results of this analysis are not extremely accurate, they do serve as a motivation to 
carry the treatment further. The reasons for the discrepancies are not clear. 
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Figure 9: Waveform of the sound generated by the prime mover at a temperature difference of 
325 °C. 
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Figure 10: Spectrum of the waveform shown in Fig. 9. 
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Figure 1 1: Waveform of the sound generated by the prime mover at a temperature difference of 
368°C. 
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Figure 12: Spectrum of the waveform shown in Fig. 11. 
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At this point, we are not in a position to determine the relative contributions from nonlinear 
propagation and thermoacoustic generation to the observed waveform distortion. At the very 
minimum, such a determination will require an accurate model for dispersion in the prime mover. 
This area is receiving attention now. Beyond that, the equations which give rise to acoustic gain 
have been linearized. A numerical solution of the coupled momentum and energy equation might 
identify significant nonlinear terms in the coefficient of gain. 

One aspect of prime movers that has been ignored until now is their ability to go into self- 
oscillation for more than one mode simultaneously. In our prime mover, the temperature 
difference required for onset for the fundamental mode is approximately 320 °C, while that for the 
first overtone (the second mode) is 465 °C. By introducing a disk with a hole of the right 
dimensions in it into the resonator at the right place, we are able increase the onset temperature 
difference for the fundamental without affecting that for the second mode. This allows us to excite 
both modes simultaneously. If the attenuator (disk) is not used for the fundamental, we are unable 
to supply enough power to the prime mover to reach onset for the second mode. The fundamental 
transports too much heat! As one might imagine, the spectrum of the waveforms generated is quite 
rich, containing harmonics of both modes as well as sum and difference frequencies. One reason 
for exciting both modes is that it gives a measure of the dispersion in the prime mover, a quantity 
needed to understand the harmonic generation discussed above. It also provides another test bed 
for theory. We intend to pursue these measurements further. 

4. Non-perfect Standine Wave Heat Engines 

It is of interest to consider the performance of thermoacoustic heat engines that operate with 
an acoustic sound field that is not a perfect standing wave. In such a system the phase between 
acoustic pressure and acoustic velocity is somewhere between zero and 90 degrees. This is a 
difficult theoretical problem to consider in any generality, given the complexity of the problem. 
The complexity can be seen in the large parameter space required to describe any specific engine. 
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There are gas parameters: sound speed, Prandtl number, Cp/c v , and mean pressure. There are 
stack geometry parameters: plate separation, plate thickness, and stack length. There are operating 
conditions: acoustic pressure amplitude, velocity amplitude, their relative phase, frequency, and 
the heat load on a heat pumping stack or the work load on a prime mover stack. These parameters 
could be expressed differently and equivalently. 

The differential equations of N. Rott 15 * 16 are capable of describing the system with all the 
detail that we need. The only significant limitation of these equations is that they are linear and are 
expected to compare well with experimental results only at low acoustic amplitudes. In some 
ways, these equations require that too many parameters and initial conditions be specified. Some 
of the parameters can be eliminated by realizing that not all of them are independent, and that the 
parameter space can be reduced to a minimal independent set by using fairly obvious 
normalizations. Two examples are: expressing the distances in the wave propagation direction in 
terms of wavelength units, and expressing distances in the orthogonal diffusion direction in terms 
of penetration depth units. However, the reduced set of parameters is still complex, and this 
interferes with the study of generalized questions such as the significance of acoustic phasing. 

In an effort to simplify the equations of Rott, Swift 2 makes additional approximations. The 
first approximation is that the length of the stack is infinitesimally short, and the second is that the 
spacing between plates in the stack is at least several boundary layers in size, such that the 
diffusive part of the thermal wave is nearly attenuated to zero at the midpoint between two plates. 

The approximation of short stack length is a reasonable way of simplifying the system, at 
least for our present purposes. For the purpose of analyzing many real or potential experimental 
systems, the approximation is not very accurate. However, for this study, there is not any 
substantial loss of information. For example, a finite length stack could be analyzed accurately by 
modeling it as many infinitesimal stacks placed in series. In fact, this is roughly equivalent to the 
method used by numerical algorithms for solving ordinary differential equations. 

The boundary layer approximation is a severe limitation for this study. The parameter space 
of interest lies somewhere between two extremes. One extreme is the "classic" thermoacoustic heat 
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engine as described by Wheatley 1 and Swift 2 . This engine operates in a sound field that has 
acoustic pressure and acoustic velocity in quadrature, and the spacing between plates of the stack is 
at least several penetration depths. The other extreme is the Stirling cycle engine, which has the 
dynamic pressure and velocity in phase, and the spacing between plates of the "regenerator" much 
less than a penetration depth. Clearly, the boundary layer approximation is not suitable for an 
analysis of engines that approach the Stirling cycle extreme. 

Our approach will be to use Swift's 17 equations initially, and then use the Rott 15 > 16 equations 
with the short stack approximation for the actual study of acoustic phasing. Swift's equations for 
acoustically stimulated enthalpy flow H2 and acoustic work flow W2 are. 
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The initial use of Swift’s equations will assist in developing the analysis structure with a simplified 
set of equations before moving to the more complicated Rott equations. Also the initial analysis 
can be compared to some of Swift's own results for purposes of error checking, and the final 
analysis can again be compared to the initial analysis in limiting cases. 

Another technique that can be applied to our problem is the graphical display of engine data, 
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using more than one independent parameter, such as a contour plot of two independent parameters. 
Figures 13 and 14 are two such contour plots based on the above equations and a purely reactive 
standing wave sound field. Figure 1 3 shows contours for a prime mover engine where the solid 
contours represent efficiency relative to Carnot (perfect) performance and the dashed contours 
represent dimensionless power density in the form of acoustic work, for a stack having a fixed 
temperature difference. Figure 14 shows contours for a heat pump where the solid contours 
represent COP relative to Camot performance and the dashed contours represent dimensionless 
power density in the form of heat flow. The parameters chosen for these plots are typical for 
engines utilizing a 10 bar helium working fluid with modest acoustic amplitudes. The two 
independent parameters are kx, which is the normalized position of the stack in the standing wave 
and T, which is the normalized temperature gradient in the stack. 

A number of interesting observations can be made from the contour plots. Most engine 
designs would not choose an operating point that coincides with the region of highest power 
density (since the efficiency is zero there), nor would one choose an operating point that coincides 
with the region of highest efficiency (since the power density is quite low there). Instead, a more 
likely choice of operating point is one that lies somewhere on the curve of best compromise 
between power density and efficiency. This curve (not shown in the figures) is the one that 
indicates the path in the T - kx plane such that a given sacrifice in efficiency, gains the most in 
jx>wer density (or vice-versa). 

The work that remains, is to produce similar contours for engine performance that are based 
on the Rott equations directly, and then vary the acoustic phasing. Also, numerically producing 
the curve of best compromise on the contour plots would be useful. 
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Prime Mover Power Sc Efficiciency Contours (3) 




Figure 13: Contour plot for a prime mover. The solid contours represent efficiency relative to 
Carnot (perfect) performance and the dashed contours represent dimensionless power 
density in the form of acoustic work, for a stack having a fixed temperature difference. 



Refrig. COPR Sc Cooling Power Contours (4) 




Figure 14: Contour plot for a heat pump. The solid contours represent COP relative to Carnot 
performance and the dashed contours represent dimensionless power density in the 
form of heat flow. 
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